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ABSTRACT 



We investigate the correlations among stellar mass (M*), disk scale length (Rd), and 



rotation velocity at 2.2 disk scale lengths (V2.2) for a sample of 81 disk-dominated galax- 



ies (disk/total > 0.9) selected from the Sloan Digital Sky Survey (SDSS). We measure 
V2.2 from long-slit Ha rotation curves and infer M* from galaxy i-band luminosities (Li) 
and g — r colors. We find logarithmic slopes of 2.60±0.13 and 3.05±0.12 for the (forward 
fit) Li — V2.2 and M* — V2.2 relations, somewhat shallower than most previous studies, 
with intrinsic scatter of 0.13 dex and 0.16 dex, respectively. Our direct estimates of the 
total-to-stellar mass ratio within 2.2Rd, assuming a Kroupa (2002) IMF, yield a median 
ratio of 2.4 for M* > 1O 1O M and 4.4 for M* = 10 9 - 1O 1O M , with large scatter at a 
given M* and Rd- The typical ratio of the rotation speed predicted for the stellar disk 
alone to the observed rotation speed at 2.2Rd is ~ 0.65. 

The distribution of scale lengths at fixed is broad, but we find no correlation 



between disk size and the residual from the M* — V2.2 relation, implying that the M* — 
V2.2 relation is an approximately edge-on view of the disk galaxy fundamental plane. 
Independent of the assumed IMF, this result implies that stellar disks do not, on average, 
dominate the mass within 2.2Rd- We discuss our results in the context of infall models 
^ ■ where disks form in adiabatically contracted cold dark matter halos. A model with 

a disk-to-halo mass ratio = 0.05 provides a reasonable match to the Rd — 
distribution for spin parameters A ranging from ~ 0.04 — 0.08, and it yields a reasonable 
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match to the mean M* — V2.2 relation. A model with = 0.1 predicts overly strong 
correlations between disk size and M* — V2.2 residual. Explaining the wide range of 
halo-to-disk mass ratios within 2.2R<i requires significant scatter in m& values, with 
systematically lower for galaxies with lower M* or lower stellar surface density £*. 

Subject headings: galaxies: photometry, kinematics and dynamics 

1. Introduction 

The tight correlation between luminosity and rotation speed is one of the fundamental charac- 
teristics of the disk galaxy population (Tully & Fisher 1977, hereafter TF). With stellar population 
modeling and HI gas measurements, this correlation can be expressed in terms of stellar mass or 
total baryonic mass, in place of luminosity (Bell & de Jong 2001; McGaugh et al. 2000). The 
form and tightness of the TF relation are critical tests for theoretical models of galaxy formation 
(e.g., Cole & Kaiser 1989; Kauffmann, White, & Guiderdoni 1993; Cole et al. 1994; Eisenstein 
& Loeb 1996; Steinmetz & Navarro 1999; Avila-Reese, Firmani, & Hernandenz 1998; Firmani <fc 
Avila- Reese 2000) . In classic models of disk galaxy formation by dissipative gravitational collapse 
(e.g., Fall k Efstathiou 1980; Gunn 1983; Dalcanton, Spergel, & Summers 1997; Mo, Mao, & White 
1998, hereafter MMW), the quantities that determine the disk rotation curve are the concentration 
parameter of the dark matter halo, the ratio of the disk baryonic mass to the total mass of the halo, 
and the disk scale length, which is determined by its angular momentum. These theoretical models 
suggest that disk size could be an important additional parameter in disk galaxy correlations (Shen, 
Mo, & Shu 2002; Dutton et al. 2005). A strong correlation between disk size and TF residual is 
also expected if disks are "maximal" and therefore make a dominant contribution to the observed 
rotation speed (Courteau & Rix 1999). 

In this paper, we investigate the correlations among rotation speed, stellar mass, and scale 
length in a sample of disk galaxies selected from the Sloan Digital Sky Survey (SDSS; York et 
al. 2000; Abazajian et al. 2004). We have obtained long-slit optical spectra of these galaxies and 
used them to extract Ha rotation curves. We use the g-r color of each galaxy to estimate its 
stellar mass-to- light ratio M*/L, and therefore its stellar mass, following the prescription of Bell 
et al. (2003). We apply bulge-disk decomposition to the SDSS i-band images to select a sample of 
disk-dominated galaxies (disk/total > 0.9) and to measure the disk exponential scale length. We 
estimate circular velocities at 2.2 disk scale lengths, where the rotation curve of a self-gravitating 
exponential disk reaches its maximum (Freeman 1970). We use the estimated stellar masses to 
separate the disk and halo contributions to the total mass within this radius, and we investigate 
the halo-to-disk ratio as a function of stellar mass and disk scale length. We discuss our results in 
the context of MMW-style disk galaxy models. 
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2. Photometric and Spectroscopic Observations 

The SDSS galaxy redshift survey has an unprecedented combination of large area, depth, and 
photometric quality, thanks to the combination of a large format camera (Gunn et al. 1998), high 
throughput multi-object spectrographs, careful calibration procedures (Fukugita et al. 1996; Hogg 
et al. 2001; Smith et al. 2002), and an efficient series of data reduction and targeting pipelines 
(Lupton et al. 2001; Stoughton et al. 2002; Strauss et al. 2002; Blanton et al. 2003a; Pier et al. 
2003; Ivezic et al. 2004). We have selected a sample of ~ 200 galaxies from the main galaxy sample 
(Strauss et al. 2002) of the SDSS redshift survey for follow-up dynamical study with long-slit Ha 
spectroscopy. Our full sample covers a representative selection of galaxies in the absolute magnitude 
range — 18 > M r > —22, with no morphological pre-selection other than an i-band isophotal 
axis ratio cut of b/a < 0.6, needed to allow accurate inclination corrections to observed rotation 
velocities. We impose a minimum redshift, cz > 5000 kms -1 , so that peculiar velocities do not cause 
large uncertainties in distance (and thus luminosity and size). We adopt a luminosity-dependent 
maximum redshift of 9000 kms" 1 (-18 > M r > -19.5), 11000 kms' 1 (-19.5 > M r > -20.5), and 
15000 km s _1 (—20.5 > M r ), so that galaxies are spatially resolved and the distribution of absolute 
magnitudes is roughly flat over the range —18 to —22. Analysis of the Tully-Fisher relation and its 
residuals for this broadly representative galaxy sample, including detailed discussion of the sample 
definition, spectroscopic observations and data reduction procedures, and rotation curve fits, will 
be presented by Pizagno et al. (in preparation, hereafter P05). Here we summarize the relevant 
aspects of these procedures and describe the selection of the disk-dominated galaxy sample that is 
analyzed in this paper. 

In a series of observing runs between June 2001 and April 2004, we obtained long-slit spectra 
covering the Ha wavelength region for a total of 234 galaxies in the velocity and absolute magnitude 
ranges described above, using the TWIN spectrograph on the Calar Alto 3.5-m telescope (189 
galaxies) and the CCDS spectrograph on the MDM 2.4-m telescope (45 galaxies). Typical exposure 
times were 30 minutes at Calar Alto and 60 minutes at MDM, with instrumental setups yielding 
FWHM resolution ~ 1.48A and ~ 1.93A respectively. We obtained usable Ha rotation curves 
for 170 galaxies, ~ 70% of the input sample, with the remainder showing insufficient extended 
Ha emission or (in a few cases) excessively irregular velocity profiles. We applied the bulge-disk 
decomposition program GALFIT (Peng et al. 2002) to the i-band images of these 170 galaxies, 
taken from the SDSS corrected frames. Specifically, we fit each galaxy with a combination of an 
inclined exponential disk and a bulge with a surface brightness profile exp[— {r /r s ) l l n ] (Sersic 1968), 
with the index n constrained to the range 0.5 > n > 5.0. 

For this paper, we select those galaxies with disk-to-total luminosity ratio fd > 0.9. These 
systems may not be perfectly described by smooth exponential disks, but the addition of a bulge 
containing more than 10% of the light does not allow a better fit. We rejected five galaxies for 
which the discrepancy between the GALFIT exponential disk position angle and the SDSS isophotal 
position angle, used for the long-slit observations, would lead to a velocity difference of more than 
10%. We confirmed the disk-dominated nature of the remaining 81 galaxies by visual inspection. 
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We note that fd > 0.9 is a stronger morphological cut than that in most TF samples, which also 
include some galaxies with significant bulges. We adopt the more stringent cut mainly because 
it allows us to define scale length and velocity measures, Rd and V2.2, that are insensitive to 
ambiguities of bulge-disk decomposition. Also, while the bulge formation mechanisms in late-type 
galaxies are uncertain, our sample gives theoretical modelers a clear target to make predictions for: 
nearly bulgeless galaxies, with an absolute magnitude distribution that is approximately flat in the 
range — 18 < M r < —22. We compare results from this sample to those of our full, morphologically 
representative sample in P05. 

We compute galaxy luminosities using SDSS Petrosian fluxes and colors using SDSS model 
colors, both K-corrected to redshift z = using Blanton et al.'s (2003b) kcorrect_v3 . lb. We 
compute distances using the SDSS heliocentric redshifts corrected to the rest frame of the Local 
Group barycenter (Willick et al. 1997), assuming a cosmological model with f2 m =0.3, $1^=0.7, and 
Hubble constant Hq = 70 km s^ 1 Mpc^ 1 . We incorporate a distance uncertainty corresponding 
to 300 km s -1 when calculating disk scale length and luminosity uncertainties, to account for the 
typical amplitude of small scale peculiar velocities (Strauss Sz Willick 1995). Figure 1 shows the 
distribution of our sample galaxies in the color-magnitude plane. Crosses show galaxies that did 
not have enough Ha emission for extended rotation curves. Although most of these failed galaxies 
are either red or low luminosity, the galaxies with successful rotation curve measurements cover all 
populated areas of the color-magnitude plane, so there are no major categories of galaxies in this 
absolute magnitude range that are missing from our sample. Filled circles show the disk-dominated 
subset analyzed in this paper. These again span all populated regions of the color-magnitude plane, 
though the distribution is somewhat bluer than that of the full sample, and the fraction of disk- 
dominated systems is higher at low luminosity. The surface brightness distribution (not shown) 
is similar to that of the full sample, though the low and intermediate luminosity galaxies that are 
eliminated by the fd > 0.9 cut tend to be above the median surface brightness. The GALFIT total 
magnitudes in i-band are 0.12 magnitudes brighter than SDSS Petrosian magnitudes on average, 
with a standard deviation of 0.10 magnitudes, in reasonable agreement with expectations (see 
Strauss et al. 2002; Graham et al. 2005). We use the Petrosian magnitudes for our analysis so that 
our results refer to quantities easily accessible from the SDSS database. 

The spectroscopic data were dark subtracted, flat-fielded, and linearized using standard IRAF 
procedures, as outlined by Massey, Valdes, & Barnes (1992). Following Courteau (1997), we extract 
2-D spectra along the spatial direction and measure the Ha emission line centroid at each location 
along the slit. The Ha emission line centroid uncertainty is between 2 and 12 kms -1 depending on 
the signal-to-noise ratio of the emission line. We define the rotation curve as the spatial variation of 
the emission line centroids along the major axis of the galaxy. We fit galaxy rotation curves with an 
arc-tangent function, which has a minimal number of parameters while still describing the global 
shape of typical rotation curves adequately (Courteau 1997). Specifically, we use a Levenberg- 
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Marquardt \ 2 minimization routine (Press et. al. 1992) to fit the data with the functional form 



where Vq is the central velocity, V c \ rc is the asymptotic circular velocity, r is the position along the 
slit, ro is the center of the rotation curve (where V = Vq), and r t is the turnover radius at which 
the rotation curve begins to flatten. When performing the fit, we add 10 kms -1 in quadrature to 
the observational error on each Ha data point, to account for non-circular motions and to ensure 
that parameters are determined by the overall shape of the rotation curve rather than the high 
signal-to-noise data points in the inner parts of the rotation curve with small uncertainties. 

We adopt the circular velocity at 2.2 disk scale lengths as our measure of rotation speed. The 
rotation curve of an isolated exponential disk peaks at this radius (Freeman 1970), and Courteau 
(1997) shows that this velocity measure produces the tightest TF relation. We infer the observed 
rotation velocity (V^if) by evaluating equation (1) at r = 2.2 x Rd, where Rd is the i-band disk 
scale length determined by GALFIT. The uncertainty in is determined using the covariance of 
the parameter errors returned by the Levenberg-Marquardt routine. We correct for inclination 
by using the GALFIT-determined disk axis ratio and the equation 



where 0.19 is the assumed intrinsic axis ratio for an edge-on disk and b/a is the measured i-band 
axis ratio. Observational estimates of the intrinsic axis ratio vary from 0.10 to 0.20 depending on 
galaxy type (Haynes & Giovanelli 1984). We choose 0.19, typical for spiral galaxies, and note that 
the range 0.10 — 0.20 corresponds to a small variation (typically ~1 kms -1 ) in V2.2- As discussed in 
detail by P05, roughly 1/3 of our galaxies have rotation curves that are still rising at the outermost 
Ha point. We have included these galaxies in our sample, but we have checked that excluding them 
makes minimal difference to our results. 

Table 1 lists the SDSS identifier, distance, i-band luminosity, g — r color, stellar mass, GALFIT 
disk fraction, disk exponential scale length, and rotation velocity V2.2 for the 81 galaxies that com- 
prise our disk-dominated sample. The luminosities and colors are corrected for internal extinction, 
and stellar masses are computed from these extinction-corrected quantities, as described in the 
next section. 



Figure 2a show the i-band TF relation for our sample of 81 disk-dominated galaxies. Following 
standard practice, we correct luminosities, and colors, for internal extinction based on the disk 
axis ratio and luminosity, using the prescription of Tully et al. (1998) interpolated to the central 
wavelength of the SDSS i-band. We convert luminosities to solar units using Mj jQ = 4.56 (Bell 




(1) 




(2) 
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Results 
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et al. 2003). The three representative error crosses in the lower right show the 90th-percentile, 
50th-percentile, and lOth-percentile values of the observational uncertainties. We add 1/3 of the 
inclination correction in quadrature to the luminosity uncertainty to represent the uncertainty in the 
inclination correction itself. Luminosity uncertainties are dominated by this inclination correction 
uncertainty and by the 300 km s -1 peculiar velocity uncertainty. 

The solid line shows our "forward" fit to the observed TF relation. Specifically, we fit a relation 

y = a(x-x )+b (3) 

with y = log Li/ L & , x = log V^/kms -1 , assuming a Gaussian intrinsic scatter of dispersion a in y 
at fixed x, in addition to observational uncertainties. We determine maximum likelihood values of 
a, b, and a by maximizing 

ln ( L ) = -\Yl ln + <v + ° 2 ^) " E 2( ^+^l V f a 2 j + constant , (4) 

where and <jj i2/ are the observational uncertainties for data point i (see P05 for discussion). We 
determine the la errors on a, b, and a by repeating this procedure for 100 bootstrap subsamples 
of the full data set and taking the dispersion as the uncertainty. We choose the reference value xo 
so that there is essentially no covariance between the errors in the slope a and intercept b. 

For the forward TF relation, we obtain 

/ y \ 2.60±0.13 

(V1O>°L ) = 1.84(±0.09) (^^) • (5) 

The intrinsic scatter is a=0.13 dex = 0.33 magnitudes, comparable to that of other TF studies 
(e.g., Courteau 1997; Kannappan, Fabricant, & Franx 2002). The dotted line shows the inverse 
TF fit, in which we assume that there is Gaussian intrinsic scatter of logV2.2 at fixed logLj instead 
of the reverse. Forward and inverse fits correspond to different hypotheses about the intrinsic 
distribution of the correlated quantities, and they yield different slopes except in the limit of zero 
intrinsic scatter. Our inclusion of the intrinsic scatter as a fit parameter means that points with 
small observational errors do not get inappropriately large weights in determining the slope and 
normalization, a difference from many previous analyses. However, for this sample the intrinsic 
scatter is small enough that the derived slope is not highly sensitive to the fitting procedure (for 
example, the inverse fit is V2.2 oc L^ 2 ' 9 ). Our slope is shallower than that found by some previous 
studies, e.g., Verheijen (2001), who finds a slope of « 4.5 in if'-band, or Kannappan, Fabricant, & 
Franx (2002), who find a slope of w 4.0 in .R-band. However, it agrees well with the r-band slope 
of 2.54 found by Courteau (1997), whose sample criteria and analysis procedures are most similar 
to ours. We discuss possible contributions to slope differences in more detail below. Our estimates 
of the slopes, intercepts, and intrinsic scatter of these and all other bivariate relations fit in this 
paper are listed in Table 2. Residuals from these relations show no discernible correlation with axis 
ratio, which indicates that our inclination corrections are accurate in the mean, even if they are not 
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perfect on a galaxy-by-galaxy basis. The Appendix presents a Monte Carlo test for Malmquist-type 
biases in our sample selection and analysis and shows that they are small compared to our quoted 
statstical errors, with an effect ~ 0.01 on the forward TF slope. 

The point types in Figure 2a encode galaxy color, relative to the expectation for the galaxy's 
i-band luminosity. We fit a linear mean relation to the sample's (g-r) vs. logLj relation and divide 
the sample into three nearly equal parts based on the residual from this relation. Squares represent 
galaxies redder than the mean by 0.02 mag or more, circles represent galaxies bluer by 0.05 mag 
or more, and triangles show the remaining galaxies. There is a slight tendency for red galaxies to 
lie below the mean TF relation, as one might expect given the higher mass-to-light ratios or red 
stellar populations, but the trend is weak relative to the scatter. 

In Figure 2b, we have converted galaxy i-band luminosities to stellar masses, using the pre- 
scription of Bell et al. (2003) to infer each galaxy's stellar mass-to-light ratio from its (g-r) color. 
We use color in preference to the spectroscopic methods of Kauffman et al. (2003) because the 
SDSS fibers cover only the central regions of these relatively nearby galaxies and may not sample 
a representative stellar population. Bell et al. (2003) adopt a "diet Salpeter" initial mass function 
(IMF) chosen so that stellar disks have the maximum mass allowed by rotation curve constraints. 
At each g-r color, we multiply their stellar mass-to-light ratios by 0.71 to correspond to a Kroupa 
IMF, which better represents direct observational estimates of the IMF (see Bell et al. 2003 for 
further discussion). Specifically, we calculate stellar masses using the relations 



where the two coefficients are from Table 7 of Bell et al. (2003). We use inclination corrected 
luminosities and colors, again based on the Tully et al. (1998) prescriptions, but the inclination 
effects are small (typically less than 10%) because galaxies move along a locus of roughly constant 
M* as their luminosities and colors are corrected for extinction. We add an error contribution to 
M* that is one-third of the applied inclination correction. 

Points in Figure 2b are again coded by galaxy g-r color residual, now computed as a function of 
stellar mass. The conversion to stellar mass has removed the slight trend of TF residual with color 
residual, as one might expect if stellar mass is the more fundamental quantity. Solid and dotted 
lines show our best-fit forward and inverse relations, which are listed in Table 1. The forward 
relation is 



The best-fit intrinsic scatter is 0.16 dex, slightly higher than 0.13 dex found for the Lj — V2.2 
relation. Some of this increase could reflect galaxy-to-galaxy variations in stellar populations or 
extinction properties, which would change the true M*/L ratios at fixed g — r color; we have 
assumed a deterministic relation between g — r and M*/L and have not included any scatter about 



M*/M = (Li/L Q ) (M*/Lj), 
log (M, /Li) = -0.222 + 0.864 x (g - r) + log 0.71 



(6a) 
(6b) 




3.05±0.12 



(7) 
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this relation in our observational error budget. The BD01 models have 0.1 dex scatter in stellar 
mass-to- light ratio (M*/Lj) at fixed color. 

The short-dashed line shows the best-fit stellar-mass TF relation found by Bell & de Jong (2001, 
hereafter BD01), using Verheijen's (2001) data for the Ursa Major cluster. We have multiplied the 
normalization of their relation (the I-band fit with mass-dependent inclination correction from 
their Table 2) by 0.71 to adjust to the Kroupa IMF assumed here. Distances are calibrated to 
Hq « 70kms _1 Mpc~ 1 in both cases, though there could be some uncertainty in the relative 
distance normalization from the peculiar velocity of Ursa Major. The two relations agree at V2.2 ~ 
200 km s -1 , but the BD01 relation is substantially steeper than ours, with a slope of 4.49 ± 0.23 vs. 
3.05 ± 0.12. At V2.2 ~ 100 km s^ 1 the BD01 relation traces the lower envelope of our data points 
and lies ~ 0.4 dex below our best-fit relations. Differences in the samples and analysis methods 
include: our use of the updated Bell et al. (2003) stellar population models in place of the BD01 
models, our use of a disk/total > 0.9 cut vs. BDOl's more generic "late type galaxy" selection, our 
use of V2.2 as a velocity measure in place of Verheijen's (2001) Va at measure used by BD01, and 
our "field" (or, more accurately, random) environment selection vs. their cluster sample, 

Since Verheijen (2001) finds a steep if '-band TF relation for the Ursa Major galaxies, BDOl's 
steep slope (relative to ours) appears to be intrinsic to the sample, not a consequence of any 
differences in stellar population modeling. Our disk/total cut also seems unlikely to be the main 
source of difference, since the I0W-V2.2 galaxies that are in our sample lie significantly above any of 
the BD01 data points, and the relatively small number of I0W-V2.2 galaxies excluded by our cut lie 
on or above the best-fit relation. To investigate the importance of velocity definition differences, 
we used the data of Courteau (1997), who lists V2.2 and the maximum velocity U max derived from 
a 5-parameter fit to optical rotation curves. The typical ratio of V max to V2.2 is higher for less 
massive galaxies, and since Vfl a t is likely to track V max more closely than V2.2, the trend goes in 
the right direction to explain the discrepancy. However, if we scale up our V2.2 values using a 
mean correction derived from the Courteau (1997) data, then our M* — V2.2 slope changes to 3.47, 
still much shallower than BD01, and the gap between the relations is still ~ 0.2 dex in M* at 
U2.2 ~ 100 km s -1 . Thus, it appears that velocity definition differences can account for roughly 
half of the difference between our results and BDOl's. We reach a similar conclusion by comparing 
Verheijen's (2001) r-band TF data points to our own (see P05). The most plausible source we can 
identify for the remaining gap is a systematic difference in properties of field and cluster spirals 
at low luminosity. Reconciling the two measurements requires the cluster galaxies to rotate faster 
by ~ 0.1 — 0.2 dex at fixed M*. Fully addressing this difference requires a large sample with a 
range of environments and both HI and optical data, so that one can mimic selection and analysis 
procedures used by different authors. 

Our stellar mass TF relation is also shallower than the baryonic TF relation derived by Mc- 
Gaugh et al. (2000), M bar /1O 1O M = 2.12(U/155.6 km s" 1 ) 3 - 98 (for H = TOkms -1 Mpc -1 ). Here 
Mbar is the sum of the stellar mass and the gas mass inferred from HI measurements. We cannot 
directly estimate Mt, ar for our galaxies because we do not have HI data, but Kannappan (2004) 



-9- 



reports a statistical correlation (with substantial scatter) between gas-to-stellar mass fraction and 
SDSS u — r color, log(G/S) = 1.46 — 1.06(u — r). If we apply this correction to our sample on a 
galaxy-by-galaxy basis, we obtain Mb ar /1O 1O M0 = 2.86(^/155.6 km s -1 ) 2,89 , shallower than equa- 
tion (7) because low mass galaxies have higher gas content. The average mass increase is 0.2 dex at 
V2.2 ~ 100 km s^ 1 and 0.05 dex at V2.2 ~ 200kms _1 . Our slope is substantially shallower than that 
found by McGaugh et al. (2000), who combine several data sets obtained in different bands. 1 The 
McGaugh et al. (2000) sample covers a much wider mass range than ours, extending to circular 
velocities V ~ 30kms _1 . The difference between optical linewidth 2V2.2 and the HI linewidth 
W20 used by McGaugh et al. (2000) could become more important at these low velocities, partly 
explaining the difference in slope, but we do not see an easy way to fully reconcile the results. 

Figure 3 plots the GALFIT disk scale length Rd against stellar mass. The dotted line shows 
the best-fit mean relation 



However, scatter about this mean relation is very broad. Points are coded by the residual from 
this best-fit relation, with squares, triangles, and circles representing the largest, intermediate, and 
smallest 1 /3 of the galaxies at a given luminosity. 

Theoretical expectations for the distribution of disk galaxies in the space of stellar mass, 
scale length, and circular velocity are clearly described by, e.g., Fall & Efstathiou (1980), MMW, 
Dalcanton, Spergel, & Summers (1997), Mo & Mao (2000), Shen, Mo, & Shu (2002), and Courteau 
et al. (2003). In this paper, we use a modeling approach similar to that of MMW to place our 
observational results in theoretical context. The central solid line in Figure 3 shows the predicted 
Rd — M* relation for galaxies with a ratio m^=0.05 of stellar mass to total halo mass formed in an 
NFW halo (Navarro, Frenk, & White 1997) with concentration parameter c=10 and spin parameter 
A=0.06. We compute this relation using equation (28) of MMW, which includes the effects of disk 
self-gravity and adiabatic contraction of the inner regions of the halo. We set the specific angular 
momentum of the disk equal to that of the halo {jd/rrid = 1 in MMW's notation), so our quantity 
A is equivalent to their A'. While the predicted relation is steeper than our best fit, it roughly 
describes the central trend of our data points. Upper and lower solid lines show the predictions for 
A=0.08 and 0.04, respectively. The envelope of these lines encloses roughly the central 80% of the 
data points. 

The distribution of halo spin parameter in N-body simulations is approximately log-normal 
with a mean (A) « 0.04 and dispersion a\ n \ s=s 0.5 (Bullock et al. 2001). For mrf=0.05, reproducing 
the Rd-M* relation requires A values in the upper half of this distribution, so disks would either 



1 McGaugh et al. (2000) also assumed a constant stellar mass-to-light ratio in each band. Had they included a 
color dependence, they would have derived a still steeper slope, since fainter galaxies are generally bluer and therefore 
have lower M*/L. 
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have to form in the higher spin halos or have a higher specific angular momentum than the dark 
matter. Since systems with low angular momentum may be more likely to form a substantial 
bulge, and thus be omitted from our disk-dominated sample, this preferential sampling of the high 
end of the A-distribution is not implausible. However, the three solid lines in Figure 3 can also 
be produced (almost exactly) with a disk-to-halo mass fraction 777^=0.025 and A of 0.03, 0.045, 
0.06; the lower rrid shifts the predicted relations to lower M*, and lower A values are required to 
compensate. Conversely, for 777(7=0. 10, the A values that yield similar Rd — M* curves are 0.055, 
0.08, 0.11. 

Figure 4 again shows the stellar-mass TF relation, M* vs. V2.2, but points are now coded by 
their residual from the best-fit R^ — M* relation, with squares, triangles, and circles representing the 
largest, intermediate, and smallest galaxies, respectively, just as in Figure 3. There is no evident 
separation among these three sets of points, i.e., no tendency of large or small galaxies to lie above 
or below the mean M* — V2.2 relation. A plot of M* — V2.2 residual against Rd — M* residual (shown 
as inset) is simply a scatter plot. As noted by Courteau & Rix (1999) and Courteau et al. (2003), 
the lack of correlation between TF residual and disk scale length argues against the "maximal disk" 
hypothesis, in which the stellar disk provides a large fraction of the rotational support at 2.2Rd, 
since in this case V22 oc GM*/ (2.2Rd). We concur with both their observational result (indeed, our 
residual correlation appears even weaker) and with their conclusion. 

Filled squares in Figure 4 show model predictions for disks with mass fraction rrid = 0.05 and 
the spin parameter A = 0.06 that yields the central solid line of Figure 3, with total halo masses 
of M h = 4 x 1O 1O M , 1.89xlO n M Q , and 1.39xl0 12 M Q . (Following MMW, we define the halo 
mass within a virial radius whose mean interior density is 200 times the critical density.) We again 
assume an initial NFW halo concentration c = 10 and compute Rd using MMWs equation (28), 
but we compute the response of the halo to the disk using the improved adiabatic contraction 
approximation of Gnedin et al. (2004), with code kindly provided by Oleg Gnedin. With these 
parameters, the model reproduces the slope and normalization of our measured M* — V2.2 relation 
as well as the Rd-M* relation. Horizontal lines attached to these points show the effect of changing 
the spin parameter to A = 0.04 and A = 0.08, corresponding to the lower and upper lines in 
Figure 3. Larger disks have weaker self-gravity and therefore lower V2.2, while compact disks make 
a substantial contribution to V2.2 and therefore spin faster. However, the predicted difference 
between large and small disks is small enough that it could plausibly be swamped by the scatter 
seen in Figure 4. 

Filled triangles in Figure 4 show models with rrid = 0.025 and A = 0.045, which also reproduce 
the central line of Figure 3, for the same three halo masses (and c = 10). Lowering rrid reduces both 

and the disk contribution to V2.2, but the shift is not exactly parallel to the M* — V2.2 relation, 
so the model predictions for rrid = 0.025 lie below the central trend of the data, by ~ 0.1 — 0.2 
dex. However, a choice of IMF with fewer low mass stars would reduce M*/Lj ratios at fixed 
color and could bring down the data points to agree with the rrid = 0.025 predictions. Because of 
the lower disk mass fraction, variations of A that reproduce the spread in the Rd — relation 
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(A = 0.03 — 0.06) produce only small shifts in V2.2 at fixed M*. 

Filled circles show models with = 0.10 and A = 0.08. The model predictions now lie above 
the central trend of the data. Adopting a more bottom-heavy IMF could raise the data points 
and remove this discrepancy, but the = 0.10 model predicts a substantial change of V2.2 over 
the range A = 0.055 — 0.11 that reproduces the spread in the Rd-M* relation. For rrid = 0.10, 
the compact galaxies (open circles) should lie noticeably to the right of the large galaxies (open 
squares) in Figure 4, and they do not. In other words, = 0.10 disks in adiabatically contracted, 
c = 10, NFW halos are too close to maximal to be consistent with the absence of a size-TF residual 
correlation. It is difficult to put this discrepancy in fully quantitative terms because a viable model 
would have to specify what parameters other than A are varying to produce the intrinsic scatter in 
the M* — V2.2 relation. We will investigate this question in future work. 

Figure 5a repackages the information in Figures 3 and 4 by plotting the rotation velocity 
predicted for the stellar disk, 



against the observed rotation velocity V2.2- The factor 0.65 in the brackets is the fraction of the 
disk mass within 2.2Rd, and the factor of 1.32 accounts for the flattened geometry of the disk 
potential (Freeman 1970; Binney & Tremaine 1987). The (1 — fd) term represents the contribution 
of the bulge, which we assume to lie entirely within 2.2i?^; this contribution is small, since fd > 0.9 
for our sample. In Figure 5a, we use the mean value of Rd at the galaxy's M*, from the best-fit 
relation shown in Figure 3 and listed in Table 1. This plot simply tilts the stellar mass TF relation 
to account for the increase of average disk size with stellar mass. The solid line shows the best-fit 
(forward) relation, 



The slope is just slightly steeper than the unit slope predicted for pure self-gravitating disks, but 
the typical offset is a factor of ~ 0.65. Since velocities add in quadrature, V^ ot = + V^, the 
normalization of equation (9) implies a typical ratio of V^/V* ~ 1.2 of the halo and stellar disk 
circular velocities at 2.2Rj. The intrinsic scatter in F*,2.2 at fixed V2.2 is 0.057 dex. If we apply 
Kannappan's (2004) color-based estimate of gas-to-stellar mass ratios to compute the total disk 
contribution V^,2.2 instead of 14,2.2) then the normalization of equation (10) rises slightly, to 0.71, 
and the slope changes to 1.07, a nearly constant ratio of disk mass to dark halo mass within 2.2Rd- 

In Figure 5b, we use each galaxy's actual scale length, instead of the mean scale length from the 
Rd — M* relation, when computing 14,2.2- The slope and normalization of the best-fit relation are 
virtually unchanged, but the intrinsic scatter is nearly two times larger, 0.105 dex instead of 0.058 
dex. Furthermore, the compact galaxies (circles) lie systematically above the mean relation (high 
14,2.2 at a given V2.2), and the large galaxies (squares) lie systematically below. Thus, even though 
the slope of the 14,2.2-^2.2 relation is close to unity, the value of 14,2.2 is, on a galaxy-by-galaxy basis, 




(9) 




(10) 
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a worse predictor of circular velocity than the stellar mass alone. Dashed lines in Figure 5 show 
V*,2.2 = O.85V2.2, which Sackett (1997) describes as a good approximation to the standard "maximal 
disk" hypothesis. 2 Our direct estimates with the Kroupa (2002) IMF lie below the maximal disk 
prediction, and they would continue to do so with the Kannappan (2004) gas correction. The 
increase of scatter between Figures 5a and 5b, another manifestation of the uncorrelated TF and 
disk size residuals, implies that this gap is not simply a consequence of underestimating M*/L 
ratios; dark halos must provide an important contribution to rotational support at 2.2Rd- 

Figure 6a plots the inferred ratio of the total mass within 2.2 disk scale lengths to the stellar 
mass within 2.2 disk scale lengths, as a function of stellar mass. The total mass is M^ 2 .2 + ^,2.2 
with 

M h , 2 . 2 = 2 -^(Vi 2 -V^ 2 ), (11) 

A^*,2.2 = 0.65/dM* + (1 — ffijM*, and given by equation (8). The total-to-stellar mass ratio 

has a flat trend with considerable scatter for > 10 10 Mq, with a median value of ~ 2.4. For 
< 10 10 Mq, the median ratio is higher (~ 4.4) in agreement with previous results (e.g., Persic, 
Salucci, & Stel 1996), and the scatter is larger. This increased mass ratio corresponds to the steeper 
than unit slope of the relation in Figure 5; the increase is reduced but not eliminated if 

we use the Kannappan (2004) gas correction to estimate baryonic masses instead of stellar masses. 
Points in Figure 6a are again coded by residual from the Rd-M* relation, and the separation of 
circles (compact galaxies) and squares (diffuse galaxies) shows that the larger galaxies, at fixed M*, 
have higher halo-to-stellar mass ratios within 2.2Rd, as one might expect. 

While the scatter in disk sizes explains some of the scatter in Mf l 2.2/^*,2.2, Figure 6b shows 
that much of the scatter must arise from another source. Here we plot 2. 2/-^*, 2. 2 against 
Rd/Rd(M*), where Rd(M*) is the mean disk scale length at the galaxy's stellar mass based on the 
best-fit Rd-M* relation (dotted line of Figure 3). Points are now coded by galaxy stellar mass, 
with M* < 1O 1O M shown by pentagons, M* > 1O 1O ' 7 M0 shown by stars, and intermediate mass 
objects shown as crosses. At a fixed Rd/Rd(M*), there is substantial scatter in Mh,2. 2/^,2. 2- The 
systematically higher halo-to-stellar mass ratios of low mass galaxies are even clearer here than in 
Figure 6a. 

The three solid curves in Figure 6b show the predictions of models with 771,^=0.05 and NFW 
halo concentrations of c = 5, 10, and 20. Along each sequence, the value of A determines the value 
of Rd, and we set Rd{M*) to be the model prediction for A=0.06, tracking the central solid line of 
Figure 3. For these calculations we have assumed M* = 10 10 M Q , but the curves are the same for 
any choice of M*. To the extent that there is a central trend of the data points, the model curves 
describe it reasonably well, but it appears that a range of halo concentrations cannot account for the 
large scatter in Mft^.2/^,2.2 at fixed disk size. The upper and lower dotted curves show predictions 



2 Specifically, Sackett (1997) states that in maximal disk decompositions, the disk rotation speed is typically 75- 
95% of the rotation speed at 2.2Rd, with the low end of the distribution populated by galaxies with large bulges, 
which would be absent from our sample. 
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for rrid = 0.025 and m^=0.10, respectively, with c=10 and Rd(M^) computed assuming A=0.045 
(for m^=0. 025) and A=0.08 (for m^=0.10). The envelope of these curves contains most of the data 
points, though there are a few with lower dark matter fractions. In the context of MMW-style disk 
galaxy models, where raj, c, and A are the parameters that determine disk properties, a substantial 
spread in rrid is required to explain the observed distribution of galaxies in the (M*, V2.2, Rd) 
space. In principle, the rotation curve shape can provide additional constraints on the disk mass 
fraction and its variation with galaxy properties (e.g., Persic & Salucci 1988; Persic, Salucci, & Stel 
1996). However, our rotation curves are not very well resolved spatially, and modeling rotation 
curve shapes requires specific assumptions about dark halo profiles, so we have not attempted to 
exploit these constraints here. 

Zavala et al. (2003) analyze the ratio of total mass to baryonic mass in a literature sample 
of disk galaxies and conclude that it correlates more directly with surface mass density than with 
luminosity or scale length individually. Figure 7 plots this correlation for our data set, using the 
mean stellar surface density within 2.2 scale lengths as the surface density measure. There is indeed 
a steady correlation over nearly two orders of magnitude in surface density, somewhat tighter than 
the correlations in Figure 6, though still with significant scatter and occasional large outliers. Most 
significantly, one can see that the total-to-stellar mass ratios for low surface density, high mass 
galaxies are similar to those of typical low mass galaxies with similar surface density. This result 
suggests, as argued by Zavala et al. (2003), that the dependences of the mass ratio on galaxy mass 
and scale length can be understood as largely reflecting a more fundamental dependence on surface 
density. The distribution of our data points in Figure 7 is similar to the distribution found by 
Zavala et al. (2003), though their sample is constructed to include more low surface brightness 
galaxies. 

Large filled squares show model predictions for rrid = 0.05, halo masses = 4 x 1O 1O M0, 
1.89x10 11 Mq, and 1.39xl0 12 M©, and spin parameter A = 0.06. Diagonal lines attached to these 
points span the range A = 0.04 — 0.08, with low-spin disks having high and low M/^2.2/-^*, 2.2- 
Filled triangles and circles show corresponding predictions for rrid = 0.025 and 0.05, with halo 
mass 1.89xlO 11 M and the A values that reproduce the three Rd — M* lines in Figure 3. The 
data roughly follow the trend predicted for A variations or nid variations as drivers of surface 
brightness variations. At fixed rrid and A, a change of halo mass does not change Mh^.2/M^ t 2.2, so 
the continuity of the trend for different mass galaxies again suggests that high mass galaxies have 
higher rrid, n °t simply higher M^. 

4. Summary and Discussion 

We have examined the correlations among stellar mass, disk scale length, and rotation velocity 
at 2.2Rd for a sample of 81 disk-dominated galaxies selected from the SDSS main galaxy redshift 
sample. The SDSS selection allows us to choose systems with a roughly flat distribution of absolute 
magnitude over the range — 18 < M r < —22 at redshifts such that peculiar velocities induce rela- 
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tively small distance uncertainties. SDSS multi-color photometry allows us to assign stellar masses 
to galaxies based on their i-band luminosities and g-r colors, using the prescription of Bell et al. 
(2003) converted to a Kroupa (2002) IMF. We use the SDSS i-band images to perform 2-d bulge- 
disk decomposition with GALFIT (Peng et al. 2002). The defining morphological characteristics of 
our sample are SDSS isophotal axis ratio b/a < 0.6 and GALFIT i-band disk-to-total luminosity 
ratios fd > 0.9. We do not apply any environmental pre-selection, so our sample should be repre- 
sentative of the range of environments in which galaxies of these morphological characteristics and 
absolute magnitudes are found. 

Our principal observational results are as follows: 

1. We find a best-fit (forward) i-band TF relation (L;/1O 1O L ) = 1.84(^.2/150 km s" 1 ) 2 - 60 , 
with an estimated intrinsic scatter of 0.13 dex, or 0.33 magnitudes. The slope is shallower than 
that found by some previous analyses (e.g., Verheijen 2001; Kannappan, Fabricant, & Franx 2002) 
but is similar to that of Courteau (1997), whose sample selection and analysis methods are closest 
to those here. Possible sources of difference include the environmental properties of the sample 
("field" vs. "cluster"), the morphological criteria (disk/total > 0.9 vs. more general "disk galaxy" 
selection), the adopted velocity measure (V2.2 vs. Vfj at or HI line width), and the fitting procedures. 
The intrinsic scatter is similar to that found for previous samples, though it rises substantially if we 
do not restrict the sample to disk-dominated galaxies (see P05). There is a weak trend for galaxies 
with redder than average colors to lie below the best-fit TF relation (low Lj). 

2. We find a best-fit (forward) stellar mass TF relation (M*/10 10 M*) = 2.32(y 2 .2/156kms- 1 ) 3 - 05 , 
with an estimated intrinsic scatter of 0.16 dex. The use of stellar mass removes the weak trend 
with color residual, suggesting that stellar mass plays a more fundamental role than luminosity in 
TF correlations. The slope is shallower than that found by BD01 for Ursa Major spirals, with good 
agreement at V2.2 ~ 200kms _1 but higher M* in our sample at V2.2 ~ 100 km s^ 1 . Possible sources 
of difference again include range of environments, morphological criteria, velocity measure, and fit- 
ting procedures. Our relation is also shallower than the baryonic mass TF relation of McGaugh et 
al. (2000), whose sample is more heterogeneous and extends to lower luminosities. For our sample, 

a statistical, color-based correction for gas-to-stellar mass fractions (Kannappan 2004) makes only 
a modest difference to the TF parameters. 

3. The R d - M* distribution has a best-fit mean relation R d = 4.0(M*/2.2 x 1O 1O M ) - 24 kpc, 
but the distribution is broad, with roughly a factor of three range in disk scale length at fixed M*. 

4. There is no discernible correlation between residuals of the R d — M* relation and residuals 
of the M* — V2.2 relation. At a given M*, compact galaxies do not rotate faster or slower than 
average. This result agrees with earlier analyses showing weak or negligible trends of TF residual 
with disk scale length or surface brightness (Zwaan et al. 1995; Courteau & Rix 1999; Verheijen 
2001). As emphasized by Courteau & Rix (1999), the lack of correlation between TF residual 
and disk size implies that disks cannot, in most cases, make a dominant contribution to rotation 
velocities at 2.2R d - This evidence for "sub-maximal" disks agrees with independent arguments 
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based on disk scale heights and vertical velocity dispersions (Bottema 1993, 1995; Kregel, van der 
Kruit, & Freeman 2005). 

5. Direct estimates of the stellar contribution to the rotation speed at 2.2i?^, based on 
the population synthesis mass-to-light ratios for a Kroupa (2002) IMF, yield a best-fit relation 
(^*, 2. 2/156 km s -1 ) = 0.64(V2.2/156kms -1 ) 1 ' 16 . Including estimated gas masses following Kannap- 
pan (2004) changes the intercept to 0.71 and the slope to 1.07. The low contribution from the 
stellar disk is consistent with the weak ^-trend noted above, and it agrees well with Kregel, van 
der Kruit, k, Freeman's (2005) estimate of 0.53 ± 0.04 based on disk scale heights. Courteau & Rix 
(1999) give an estimate of 0.6, also in good agreement with our results. The scatter between K,2.2 
and V2.2 is larger (0.11 dex vs. 0.06 dex) if we use each galaxy's actual Rd instead of the mean Rd 
at the galaxy's M*, another sign of the weak correlation between size and rotation speed at fixed 
M*. 

6. The ratio of halo-to-stellar mass within 1.2R& has a large range at a given M*. Median 
ratios are 2.4 for galaxies with M* > 10 10 M Q and 4.4 for 1O 9 M < M* < 10 10 M o . At any M#, 
compact galaxies have lower M/^2.2/^,2.2 and large galaxies have higher M^^.2/^,2.2, as expected. 
However, the spread in radius accounts for only a fraction of the scatter in halo-to-stellar mass ratio. 
Galaxies with a wide range of M* and Rd trace out a continuous correlation between halo-to-stellar 
mass and disk surface density, in agreement with Zavala et al. (2003), but the scatter about the 
mean trend is substantial. 

For early-type galaxies, the "fundamental plane" is close to a virial relation for the stellar 
component, a 2 ~ GM*/i? c g (Djorgovski & Davis 1987; Dressier et al. 1987; Bernardi et al. 2002), 
and the scatter of the bivariate L — a relation is significantly larger than the scatter about the 
fundamental plane. The near-virial form is naturally explained if stars dominate the central gravi- 
tational potential that determines the observed velocity dispersion, with only a modest contribution 
from dark matter (see Rusin & Kochanek 2005 and references therein). For disk galaxies, our re- 
sults show the opposite: V2.2 is better correlated with (or L) than with GM*/Rd, making the 

— V2.2 an essentially edge-on view of the disk galaxy fundamental plane. The minimal effect of 
disk size on V2.2 implies that dark matter must contribute a large fraction of the mass within the 
central two scale lengths of disk galaxies, in accord with our direct (but IMF-dependent) inference 
based on stellar mass-to-light ratios. Our conclusions on these points agree with those of Courteau 
& Rix (1999) and with the more recent analysis of Courteau et al. (2003), who investigate the 
scaling relations of disk galaxy properties in a larger but less tightly defined sample. 

To put our observational results in context, we have compared them to the predictions of theo- 
retical models in which disks form by the dissipative collapse of gas in cold dark matter halos (Fall 
& Efstathiou 1980; Dalcanton, Spergel, & Summers 1997; MMW; Shen, Mo, & Shu 2002; Dutton 
et al. 2005). The scale length and rotation speed of a galaxy with specified are determined 
in these models by the ratio of disk mass to halo virial mass (m^), the spin parameter (A), and 
the concentration parameter (c) of the NFW halo. In practice, variations of c within the expected 
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range have only moderate impact on Rd and V2.2, in part because adiabatic contraction alters the 
inner regions of the dark matter halo. A model with rrid = 0.05 reproduces the observed Rd — 
distribution if A values range from ~ 0.04 to ~ 0.08, a span that omits the lower half and the 
extreme upper tail of the log-normal A distribution predicted for dark matter halos (e.g. Bullock 
et al. 2001). For rrid = 0.025 and rrid = 0.10, the required values of A are, respectively, lower by 
25% and higher by ~ 30%. 

The rrid = 0.05 model reproduces our measured M* — V2.2 relation reasonably well, and it 
predicts a weak correlation between size and TF residual that could plausibly be washed out by the 
TF scatter. (Note, however, that such a model with fixed rrid does not explain the magnitude of 
the inferred intrinsic scatter.) The rrid = 0.025 model predicts lower M* at a given V2.2, so it would 
require a different IMF (with fewer low mass stars) to be consistent with our data. The rrid = 0.10 
model predicts slightly higher M* values at a given V2.2, and it predicts substantial size-TF residual 
correlations because of its high mass disks. These strong residual correlations appear incompatible 
with our data. 

The model comparisons in Figure 6b suggest that no model with a single value of rrid will 
reproduce our observed distribution of data points: the large scatter in M^ 2 .2/-^*,2.2 at fixed 
Rd/Rd(M*) can only be explained with scatter in the ratio of stellar mass to total halo mass. For 
our Kroupa (2002) IMF normalization, Figure 6b suggests rrid values in the range ~ 0.025 — 0.05 
for galaxies with 10 9 M o < M* < 1O 1O M , and ~ 0.05 - 0.10 for galaxies with M* > 1O 1O M . 
Systematically lower rrid values for low mass galaxies are a plausible signature of more efficient 
supernova feedback in shallow gravitational potential wells (Dekel k, Silk 1986) or in galaxies that 
lack a confining envelope of shock-heated gas Keres et al. (2004). 

We have not attempted to develop a complete and self-consistent model for the distribution 
of disk mass functions and halo parameters needed to reproduce the observed, joint distribution of 
M*, Rd, and V2.2- The results above suggest the outline of such a model: disk mass fractions span 
a substantial range, with a central value rrid ~ 0.05 and systematically lower values for low mass 
galaxies, and the disk galaxy population samples mainly the upper half of the halo spin parameter 
distribution. Our observations provide the data needed to constrain more complete models of this 
sort, along the lines pursued by Shen, Mo, & Shu (2002) and Dutton et al. (2005), or to test 
the predictions of hydrodynamic simulations of disk galaxy formation (e.g. Navarro & Steinmetz 
1997). Confrontations between these models and our data should lead to better understanding of 
the mechanisms that govern disk galaxy formation and the relations between dark and luminous 
matter in the inner regions of disk galaxies. 
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A. Appendix 

Random errors in galaxy distances or magnitudes can lead to biased estimates of TF parameters 
because there are more distant and faint galaxies to scatter in one direction than vice versa, and 
because objects can scatter across selection boundaries. Most discussions of these "Malmquist"- 
type biases have focused on apparent magnitude limited samples, or on the systematic bias in the 
derived peculiar velocity field (e.g., Lynden-Bell et al. 1988; Gould 1993; Teerikorpi 1993; Strauss 
& Willick 1995). Since our selection procedure is quite different from those of most previous TF 
surveys, we have tested for Malmquist-type biases with a simple Monte Carlo experiment. 

The measurement errors in SDSS Petrosian magnitudes are generally small, and luminosity 
uncertainties are therefore dominated by line-of-sight peculiar velocities, which we have assumed in 
our analysis to be drawn from a Gaussian of dispersion 300 km s -1 . For our Monte Carlo sample, 
we assign random 3-d positions to 10 6 artificial galaxies and draw their absolute r-band magnitudes 
from the SDSS luminosity function of Blanton et al. (2003c). We assign each of these galaxies a 
circular velocity drawn from the best-fit i-band inverse TF relation listed in Table 2, with slope 
a = 0.346, intercept b = log V^/kms -1 = 2.185 at log L/L & = 10.293, and intrinsic scatter of 
a = 0.048 dex in log Vi.i- For purposes of this experiment, we ignore the slight difference between 
i-band and r-band TF relations, since we only need a qualitatively realistic assignment for the test. 

We modify each galaxy's redshift by a peculiar velocity drawn from a 300 km s -1 Gaussian, 
and we apply the same absolute-magnitude dependent redshift cuts that we used for our sample 
definition (see §2), using the apparent rather than true absolute magnitude. Finally, we randomly 
draw from this cut sample a subset of galaxies that matches the nearly flat M r distribution of 
our observed sample, selecting ten times as many artificial galaxies as observed galaxies in each 
absolute magnitude bin. Applying our maximum likelihood estimation method to this artificial 
sample yields inverse TF parameters a = 0.347, b = 2.184, and a = 0.050, in excellent agreement 
with our input values. The forward TF parameters are a = 2.613, b = 10.26, and a = 0.14 
mag, in excellent agreement with the values of 2.603 ± 0.133, 10.266 ± 0.020, and 0.131 ± 0.015 
derived for our observed sample (see Table 2). We conclude that any Malmquist-type biases in our 
sample definition or analysis are smaller than our quoted statistical errors. This small impact is not 
surprising, since the fractional distance errors are small and the velocity and absolute magnitude 
range of the sample are fairly large, making scatter across selection boundaries a minimal effect. 
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Note. — Local Group barycenter distances are in Mpc, for Ho = 70kms~ 1 Mpc -1 . Luminosities are derived from SDSS 
Petrosian magnitudes and corrected for internal dust extinction following Tully et al. (1998). Colors are SDSS model colors 
corrected for internal extinction. The observed Petrosian magnitudes, model colors, helio-centric redshifts, and corrected frames 
can be obtained from the public SDSS DR2 server at http://www.sdss.org/dr2. 

a GALFIT magnitude used instead of SDSS Petrosian, because of poor pipeline estimation of the total observed magnitude. 
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Table 2. Bivariate Relation Fits 
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Note. — Bivariate relations are fit with the model y=a(x-xo)+b with a Gaus- 
sian intrinsic scatter of y at fixed x, with dispersion a. Errors listed for a, b, and 
a are computed from 100 bootstrap trials, and the value of xo is chosen so that 
errors in a and b are uncorrclatcd. Luminosities are in Lj,©, stellar masses in Mq, 
velocities in kms -1 , and radii in kpc. 

a V».2.2 is computed using the mean disk scale length for the galaxy's stellar 
mass. 

b V*,2.2 is computed using the each galaxy's measured scale length. 
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Fig. 1. — Distribution of sample galaxies in the color-absolute magnitude plane. Filled circles are 
the sample of 81 disk-dominated galaxies analyzed in this paper, and open circles are galaxies with 
fd less than 0.9. Crosses show galaxies with insufficient extended Ha for useful rotation curve 
measurements. 
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Fig. 2. — faj TF relation of our disk-dominated sample: inclination corrected i-band luminosity 
vs. circular velocity at 2.2 disk scale lengths. Circles, triangles, and squares show galaxies that are 
blue, intermediate, and red with respect to the mean color-absolute magnitude relation. Solid and 
dotted lines show the best forward and inverse fits to the data points (see Table 2 for parameters). 
(b) Same as (a), but with Lj replaced by the estimated stellar mass M*. Here and in subsequent 
figures, M* is estimated from the z-band luminosity and a stellar mass-to-light ratio inferred from 
the g-r color following Bell et al. (2003) (see equation 6). The short-dashed line shows the best-fit 
stellar mass TF relation from BD01, adjusted to the Kroupa (IMF). In this and all subsequent 
figures, error crosses show the 90th, 50th, and lOth-percentile values of the observational errors. 
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log M^/M 

Fig. 3. — Relation between disk scale length and stellar mass. The dotted line is the best-fit 
mean relation, and points are coded by distance from this relation. The central solid line shows 
the predicted relation for exponential disks formed in adiabatically contracted NFW halos with 
concentration c = 10, spin parameter A = 0.06, and ratio = 0.05 of disk mass to halo virial 
mass. Lower and upper solid lines show the model predictions for A = 0.04 and 0.08, respectively. 
Other combinations of (md,A) can produce similar results, as discussed in the text. 
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Fig. 4. — Stellar mass TF relation, as in Fig. 2b, but with points coded by their residual from the 
mean Rj — M* relation, as in Fig. 3. The inset panel plots TF residual against residual. In the 
main panel, filled squares show predictions for a model with md = 0.05, c = 10, and A = 0.06, for 
three different halo masses. Horizontal error bars on these points show the change in V2.2 when A is 
changed to 0.04 (higher V2.2) or 0.08 (lower V2.2)) the values corresponding roughly to the range of 
observed disk sizes. Filled circles and filled triangles show corresponding predictions for rrid = 0.10 
and 0.025, respectively, with A values chosen to produce the observed range of disk sizes for the 
corresponding (see text for further discussion). 
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Fig. 5. — Comparison of the rotation velocity predicted for each galaxy's stellar component 

(eq. [8]) to the observed Va.2- Points are coded by residual from the Rd — M* relation, as in Fig. 3. 
In panel (a), K,2.2 is computed using the mean Rd at each galaxy's M*, while in panel (b) it is 
computed using each galaxy's observed Rd- Solid and dotted lines show forward and inverse fits 
with parameters listed in Table 2. The dashed line shows the relation V* 2.2 = O.85V2.2 expected 
for "maximal" disks (Sackett 1997). 
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Fig. 6. — (a) Ratio of total mass within 2.2i?^ to stellar mass within 2.2i?^, as a function of M*. 
Points are coded by residual from the Rd — M* relation. Filled circles with error bars show the 
mean and standard deviation for galaxies in the mass ranges log M^/Mq > 10.7, 10 — 10.7, and 
< 10. (b) Similar to (a), but mass ratios are now plotted against the ratio of each galaxy's scale 
length to the mean scale length for its stellar mass. Crosses, stars, and pentagons show galaxies 
with log M*/M & > 10.7, 10 — 10.7, and < 10, respectively. The three solid curves show predictions 
of models with = 0.05 and concentration parameters c = 5, 10, and 20 (bottom to top). The 
dotted curves show predictions for rrid = 0.025 (upper) and = 0.1 (lower), for c = 10. 
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Fig. 7. — Total-to-stellar mass ratio as a function of the mean stellar surface density within 2.2Rd- 
Stars, crosses, and pentagons represent high, intermediate, and low mass galaxies, as in Fig. 6b. 
The filled triangle, square, and circle with attached diagonal lines represent models with (rrid, A) = 
(0.025,0.045), (0.05,0.06), and (0.10,0.08), respectively, with c = 10 and a halo mass M h = 1.89 x 
1O 11 M0. The points correspond to the three central models shown in Fig. 4, and lines show the 
effect of varying A over the range that reproduces the observed range of R^, as in Fig. 4. Additional 
filled squares show the m d = 0.05 model for M h = 4 x 1O 1O M (low £*) and 1.39 x 1O 12 M (high 
£*); the lines of varying A would approximately parallel those of the central model. 



